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Abstract The richness of phenomena occurring in heavy-ion colHsions calls for 
microscopic approaches where the motion of each nucleon is treated quantum me- 
chanically. The most popular microscopic approach for low-energy collisions be- 
tween atomic nuclei is the time-dependent Hartree-Fock (TDHF) theory, providing 
a mean-field dynamics of the system. The TDHF approach and some of its exten- 
sions are used to predict the evolution of out-of-equilibrium nuclear systems. The 
formation of di-nuclear systems with a structure close to molecular states is in- 
vestigated. In particular, lifetimes and exit channels are described. The formation 
of light molecules and the dynamics of a-clustering are discussed. Di-nuclear sys- 
tems formed in transfer, deep-inelastic, and quasi-fission reactions, as well as hyper- 
heavy molecules produced in reactions between actinides are also investigated. The 
formation and stability of structures in neutron star crusts are finally discussed. 



1 Introduction 

Clustering in atomic nuclei is a general concept which includes a large variety of 
phenomena. Most of them have been covered in the volumes of "Clusters in Nuclei". 
These include, for instance, a-clustering [1, 2, 3, 4, 5] and molecules formed by two 
light [6, 7, 8], intermediate/heavy [9] or very heavy [10, 3] fragments. 

These cluster configurations are usually considered as specific structures of the 
total systems. However, except for light nuclei such as some beryllium isotopes 
[1, 11,3], and eventually some heavier nuclei subject to cluster radioactivity [12, 3], 
they are barely found in nuclear ground-state but rather in (sometimes highly) ex- 
cited states. The question of the formation of these systems in nuclear reactions 
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comes then naturally. However, the large variety of structures and reactions to in- 
vestigate makes it very challenging to develop a unique model able to describe both 
dynamical and static properties of these systems. 

One possibility to describe clustering structures and dynamics is to assume a 
priori the presence of clusters in the state of the system. This is done, for instance, 
to describe a-condenstates in Ref. [5] and di-nuclear systems in Ref. [9]. Combining 
microscopic and macroscopic degrees of freedom is also possible within the two- 
center shell model [12, 10] and the generator coordinate method (GCM) [4]. 

Nevertheless, purely microscopic approaches describing the state of each nu- 
cleon can be used assuming specific forms of the many-body wave functions. In the 
anti-symmetrised molecular dynamics (AMD) model, for instance, Gaussian single- 
particle wave-functions are considered [1, 11]. Another example of purely micro- 
scopic approach is the time-dependent Hartree-Fock (TDHF) theory (see Ref. [13] 
for a recent review), which will be thoroughly used in the present chapter. 

An attracting feature of the TDHF approach is that it uses the same energy den- 
sity functional for both the structure of the collision partners and their dynamics 
during the reaction. Also, the same functional is used over the nuclear chart, al- 
lowing for both descriptions of structures and reaction mechanisms with a limiting 
number of parameters. In particular, only nuclear structure inputs are used in the 
fitting procedure of the functional. 

These aspects of the TDHF theory make it a promising tool to investigate vari- 
ous interplays between nuclear dynamics and (at least some) clustering effects. Few 
early TDHF codes have been used to study some cluster states and molecular struc- 
tures. For instance, a-clustering were investigated [14], as well as light [15] and 
hyper-heavy [16, 17] molecules. 

More recently, the dynamics of di-nuclear states formed in heavy-ion collisions 
were analysed with modern three-dimensional TDHF codes. For instance, the path 
to fusion and nucleus-nucleus potentials have been studied in Refs. [18, 19, 20, 21, 
22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 7, 
43, 44, 45, 46, 13, 47, 48]. The transfer of one or many-nucleons and the isospin 
equilibration between the fragments in contact have also been investigated in Refs. 
[19, 27, 28, 29, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 44, 61, 13, 47, 48]. 
In addition to heavy-ion collisions, the TDHF approach has been used to describe 
neutron star crust dynamics [62, 63]. 

First, formal and practical aspects of TDHF calculations are presented in Sec. 2. 
Then, the formation of light molecules and the dynamics of a-cluster states are dis- 
cussed in Sees. 3 and 4, respectively. Quasi-elastic transfer is the subject of Sec. 5, 
followed by the study of more damped collisions, namely deep inelastic collisions in 
Sec. 6 and quasi-fission reactions in Sec. 7. Calculations of hyper-heavy molecules 
dynamics in reactions between actinides are presented in Sec. 8. Finally, recent 
TDHF studies of neutron star crust dynamics are discussed in Sec. 9. 
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The TDHF theory has been developed by Dirac in 1930 [64]. This is an extension 
of the mean-field approach to the ground-state of many-fermion systems introduced 
by Hartree [65] and Fock [66]. 



2.1 The mean-field approximation 

The TDHF theory determines the dynamics of a many-fermion system out of equi- 
librium under the approximation that the state of the system can be described by 
an independent-particle state at any time. The spatial correlations between the par- 
ticles are obtained from the self-consistent mean-field. It is then assumed that each 
particle evolves independently in the mean-field generated by all the others. 

The TDHF approach is naturally well adapted to many-body systems in weak 
interactions. Indeed, when the interactions are strong, the system is expected to de- 
velop correlations which make the independent particle picture fail on a relatively 
short time scale. One may then wonder why the TDHF approach has been so suc- 
cessful in describing low-energy nuclear dynamics (see Refs. [67, 13] for reviews). 

In fact, at low energy, the Pauli principle prevents collisions between nucleons 
in such a way that the mean-free path of a nucleon in the nucleus is of the order 
of the size of the nucleus. This means that the wave-functions of the nucleons are 
essentially sensitive to the mean-field directly determined by the density. As a result, 
a nucleus in its ground state, where all the single-particle states below the Fermi 
level are almost entirely occupied, can be described with an independent particle 
state in a first approximation. Hartree-Fock calculations based on energy density 
functionals (EDF) are indeed able to reproduce quite well the binding energies and 
ground-state deformations along the nuclear chart (See Ref. [68] for a review). 

Similarly, low-energy heavy-ion collisions can be treated at the mean-field level. 
Indeed, at energies around the fusion barrier, the motion of the nuclei is slow enough 
to prevent nucleon-nucleon collisions thanks to the Pauli principle during the first 
few zeptoseconds (zs) of the reaction. However, the mean-field approximation is 
expected to fail at higher energies (e.g., in the Fermi regime), or for longer times. 
For instance, only the first steps of the fusion process can be described with TDHF, 
while beyond mean-field correlations are needed to form an equilibrated compound 
nucleus (CN) on a longer time scale. 



2.2 Formalism 

The time-dependent Hartree-Fock equation reads [64] 
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ihj^p{t)^[h[p{t%p{t)], (1) 

where p (f ) is the one-body density matrix with matrix elements 

Pafi^{^\alaa\^). (2) 

As the system is described by an independent particle-state, the state | is a Slater 
determinant of the form 

where A is the number of particles, a] creates an occupied single-particle state = 
a\\—), and | — } is the vacuum. The one -body density matrix associated to such a state 
is a projector onto the subspace of occupied single-particle states with the property 
= p. It is associated to the single-particle operator 

A 

The Hartree-Fock single-particle Hamiltonian h[p] entering the TDHF equa- 
tion (1) is obtained from the first derivative of the expectation value of the many- 
body Hamiltonian H according to the one-body density matrix. Its matrix elements 
read 

r n 8(<P\H\^) 

op^a 

In nuclear physics, however, the hard-core of the nucleon-nucleon interaction leads 
to a divergence of {^\H\'t>) when evaluated on an independent-particle state. The 
effect of the hard-core can be renormalised in such a way that the energy does not 
diverge without affecting the low-energy properties of the system. In practical HP 
and TDHF calculations, the expectation value of H is replaced by an energy den- 
sity functional E[p] determined, for instance, from the Skyrme [69] or Gogny [70] 
phenomenological effective interaction. The HF Hamiltonian then reads 

ma,= '-^. (6) 



2.3 The Skyrme energy density functional 

The Skyrme energy density functional is obtained from a zero-range effective in- 
teraction with gradient terms [69] . Numerical calculations are greatly helped by the 
zero-range nature of the interaction which simplifies the expression of the mean- 
field. 
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The interaction between two nucleons in the medium reads 

v(l,2) -ro (l+xoPa) 8 



t2 {1+X2pa) {k'-Sk 
i ?3 (1 +X3pa) P"(R) 5 

iWoa-(k'x5k) (7) 



where 5 = 5(f(l) -f(2)), k = ^(V(l)-V(2)) (relative momentum), k' is the 
complex conjugated of k acting on the left, and R = (f(l)+f(2)) /2. The op- 
erators a ~ O'(l) + d{2), with d{i) = dx{i)ex + dy{i)ey + d';(i)e;, are expressed 
in terms of the Pauli matrices ^^/^/^(O acting on the spin of the particle /. Pa = 
[1 + a(l) • d{2)] /2 corresponds to the exchange of the spins. The particle density 
in r is noted p(r) = ^j^p(ri^,ri^) where p is the one-body density matrix, s the 
spin and q the isospin. The "fi" and "f2" terms are non-local in space and simu- 
late the short range part of the interaction. Finally the last term accounts for the 
spin-orbit interaction. 

The EDF describes how the energy of the system depends on its density. It is 
standard to introduce a local energy density J^{r) such that 



E = J rfrJf(r). (8) 
This energy density can be expressed as [71] 

+B9(pV-J+j-VxS + J2p./V 



+fiioS2+Bn^S2+Bi2p«S2+Bi3p"52S^ 



q q 



(9) 



The densities entering Eq. (9) are the local density 
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where s denotes the spin, and 9, are occupied single-particle states, the kinetic en- 
ergy density 

T(r)=^ |V(p,(r.)|2, (11) 

is 

the current density 

j(r) = :^E'?''*(r^) V(p,(o) +c.c., (12) 

where c.c. means "complex conjugated", the gradient of the spin-orbit density 

V.J(r) = -;XV(p*(ri) x V(p,-(r/) • (i|c7|/), (13) 

and the spin density 

S{r)=Y^(p*irs)(piirs'){s\a\s'). (14) 

is 

In the above expressions, it is assumed that the one-body density-matrix p is diag- 
onal in isospin. The isospin is then omitted to simplify the notation. The j and S 
densities are time-odd and vanish in time-reversal invariant systems. They are, how- 
ever, important in time-dependent calculations to ensure Galilean invariance [72]. 

Note that the general Skyrme EDF includes other terms which are neglected in 
Eq. (9). These additional terms are of the form S • AS and with other densities, i.e., 

the spin-current pseudo-tensor J and the spin-kinetic energy density T [72, 22]. 

They are sometimes included in TDHF calculations [22, 21, 73]. In Eq. (9), only 

o 

the anti-symmetric part of J , which is the spin-orbit density J, is included. The 
spin-orbit energy is indeed expected to be more important (by about one order of 
magnitude) than the other spin-gradient terms [74] . 

The coefficients B, in Eq. (9) are related to the parameters of the Skyrme effective 
interactions {fo_3,;ico_3,a,Wo} as 
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(15) 



The Skyrme-HF mean-field is derived from Eq. (6). Its action on single-particle 
wave functions is then given by [71] 



{h[p](pi){r,s) = 



2'«;,(r) 



V+[/,,(r)+'C,,(r)-V 5,, 



-V,,(r).(^|(T|^') + 'W,,(r).((^|(T|^')xV) 



(16) 



where qi is the isospin of the state \^)\). The derivatives act on each term sitting on 
their right, including the wave function. The fields (functions of r) read 
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2 m* 2 m 



B^p+B^Pg (17) 



Uq = 2Bip+2B2pc,+B3{x + iV-j)+B4{z,, + N-j^) 
+2B5Ap+2B6Ap,, + {2 + a)Bjp^+" 

+B,[ap"-'Y,p^^+2p"p,]+Bg{V-J + V-J,) 

<? 

+ap«-i(Bi2S2 + Bi3l^s2) (18) 

= BgV X (j+jy)+2BioS + 2BnSq 

+2p«(Bi2S + Bi3S^) (19) 

W, = -Bgy{p+p„) (20) 

= 2B3j + 2B4j^-B9Vx(S + S,,), (21) 

where the derivatives act on the first term sitting on their right only. The label q 
denotes the isospin. The effective mass m* of nucleons with isospin q is introduced 
in [Eq. (17)]. It originates from the non-local terms of the effective interaction in 
Eq. (7). 

The parameters of the Skyrme EDF are fitted on few quantities (see, e.g., 
Ref. [74]). These include the density po ~ 0.16 fm^^ and energy per nucleon 
E /A —16 MeV of the infinite symmetric nuclear matter at saturation as well as its 
compressibility. Depending on the parametrisation, the equation of state of the infi- 
nite neutron matter [75], the enhancement factor of the Thomas-Reiche-Kuhn sum 
rule, the symmetry energy, and the radii and binding energies of few doubly-magic 
nuclei may be included as constraints into the fitting procedure as well. 

It is interesting to note that no input on nuclear reaction mechanism, such as 
fusion barriers or cross-sections, are included in the fit. Nevertheless, as we will see 
in this chapter, the description of collision dynamics with TDHF is very realistic, 
and agreements with experimental observables are sometimes impressive. 

Finally, the Coulomb interaction between the protons is added to the Skyrme 
mean-field. The direct part of the Coulomb energy reads 

,dir_e^ /■ ,3„ /■ ,3 /Pp(r)pp(r') 



Er-^Jd^rJ,VO^^. (22) 

The latter is usually computed by solving, first, the Poisson equation to get the 
Coulomb potential Vc{r), and, then, by evaluating the integral j J cPrppVc- The ex- 
change part of the Coulomb energy is usually determined within the Slater approx- 
imation as 



d'rp„{r)T. (23) 
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As a result, the contribution of the Coulomb interaction to the proton mean-field 
reads 
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2.4 Numerical implementation 



The TDHF equation is never implemented with its Liouville-von Neumann form 
given in Eq. (1). Instead, a set of non-linear Schrodinger-like equations for single- 
particle motion is used. Indeed, Eq. (1) can be expressed in a fully equivalent way 
as 

inj^Wi{t))=h\p{t)]W>{t)). 1 < i<A. (25) 

The numerical advantage of using Eqs. (25) instead of Eq. (1) is obvious in terms 
of computer memory. Indeed, for a basis of single-particle states (e.g., the number 
of points of a cartesian grid), the storage of p requires aNxN array, while the wave- 
functions in Eq. (25) require a N xA array. As A is usually much smaller than A^, 
one clearly realises the advantage of solving Eqs. (25) instead of Eq. (1). 

Equations (25) are coupled by the self-consistency of the HF Hamiltonian as 
it depends on the total density of the system. As a result, the HF Hamiltonian is 
time-dependent and one needs to solve Eqs. (25) iteratively in time. The sates at 
time t + At are determined from the states at time t assuming that h is constant 
between t and f + At. This implies that At has to be chosen small enough for this 
condition to be valid. Typical time step increments z\f in nuclear physics range from 
~ 5 X 10-25 s [76, 77] to - 1.5 X 10-^4 s [18, 78]. 

In addition, to conserve energy and particle number, the algorithm has to be sym- 
metric under time-reversal transformation. This implies that the HF Hamiltonian has 
to be evaluated at f + ^ [79]. The evolution operator then reads 



|(p,-(f + 4f))~/('- 



--'/^'l<P,(0). 



(26) 



A truncated Taylor development of the exponential is usually considered. The evo- 
lution operator then breaks unitarity and one should check the orthonormalisation 
of the wave-functions during the time evolution. 

A possible algorithm to perform the time evolution in Eq. (26) is described be- 
low: 



{I-Ph-I^f)} 

1l 



p("+i). 



4 



(27) 



p(n)+p(«+l) 
2 



p("+l) 
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where |(P; ) is an approximation of \(pi{t„ = iiAt)). A first evolution over Zif is 
performed to estimate the density at t + At. The latter is used, together with the 
density at t, to determine the density, and then the HF Hamiltonian, at f + This 
Hamiltonian is finally used to evolve the wave-function from r to f + At. 

Possible single-particle bases to solve the TDHF equation numerically are the 
harmonic-oscillator basis [80], basis-spline collocation method [81], wavelets [62], 
adaptive networks [76], or regular cartesian grids [18, 82]. Typical regular mesh 
spacing with Ax ~ 0.6 fm [77], 0.8 fm [18] and 1.0 fm [78, 82] are used. 

The initial condition of a TDHF calculation of heavy-ion collisions usually as- 
sumes that the nuclei are at some finite distance in their HF ground state. HF cal- 
culations of the collision partners then need to be performed prior to the TDHF 
evolution. This is done with the same EDF as in the TDHF calculation to ensure full 
self-consistency between structure and dynamics. Large initial distances between 
the centers-of-mass should be used to enable a proper treatment of the Coulomb 
excitation in the entrance channel. Typical distances of the order of ^ 40 fm are 
considered as a good compromise to limit computational time. It is also usually as- 
sumed that the nuclei followed a Rutherford trajectory prior to this initial condition. 
It determines the initial momenta Tik to be applied to the nucleons using Galilean 
boosts of the form 

|(p,(r = 0))=e-'«|(pf^), (28) 

where \(pf^^) are the HF single particle states. 

More details on numerical implementations of the TDHF equation can be found 
inRefs. [83, 13]. 



2.5 Beyond the TDHF approach 

The independent-particle approximation can be considered as a zeroth order ap- 
proximation to the many-body problem. In fact, the exact evolution of the one-body 
density-matrix reads 

iTi^p, = [fi,pi] + ^Tr2[vi2,pi2], (29) 

where vn is the antisymmetrised two-body interaction and pi and pi2 are the one- 
and two-body density-matrices, respectively. Solving Eq. (29) requires the knowl- 
edge of Pi2(0- The latter obeys an evolution equation which depends on the three- 
body density-matrix. In fact, Eq. (29) is the first equation of the BBGKY hierar- 
chy [84, 85, 86] providing a set of coupled equations for pi, pi2, Pi23--- 

We see that the TDHF equation is obtained by neglecting the last term in Eq. (29). 
It is important to know what is the physical meaning of this term. It contains the so- 
called two-body correlations which develop because of the residual interaction, i.e., 
the difference between the exact and mean-field Hamiltonians. Three main types of 
correlations can be identified: 
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• pairing correlations, 

• correlations induced by a collision term, 

• and long-range dynamical fluctuations. 

Pairing correlations are important for a proper description of mid-shell nuclei, 
as well as to describe pair-transfer reactions, as we will see in Sec. 5.2. They can 
be included using a "generalised" mean-field approximation. In this case, the state 
of the system is described as a quasi-particle vacuum [87]. This leads to the BCS 
model for pairing between time-reversed states, or, more generally, to the Hartree- 
Fock-Bogoliubov approximation. Nuclear dynamics in presence of pairing has been 
investigated recently with the TD-BCS approach [88, 89], and at the TDHFB level 
[90, 91, 80]. The linearised version of TDHFB is the quasi-particle random-phase 
approximation (QRPA) which has been widely applied to study nuclear vibrations 
[92, 93, 94, 95]. Applications of the TDHFB formalism to study pairing vibrations 
are presented in Sec. 5.3. 

The collision term is important at high energy, and to describe long-term dynam- 
ics such as the thermalisation of the compound nucleus. It can be added to the TDHF 
equation in what becomes the Extended-TDHF formalism [96, 97, 98, 99, 100, 101]. 
A more general approach including pairing and a collision term is given by the time- 
dependent density-matrix (TDDM) formalism [102, 103] which has been applied to 
heavy-ion collisions [104, 105]. However, to describe reaction mechanisms at inter- 
mediate energy such as multi-fragmentation, quantum effect can be neglected in a 
first approximation. Semi-classical versions of the mean-field theory including col- 
lision terms, such as the Landau-Vlasov formalism, have then been widely used to 
describe reaction mechanisms at intermediate energy [106, 107, 108]. 

Long-range dynamical fluctuations may play an important role even at low en- 
ergy. For instance, they are crucial to determine fluctuations of one-body operators, 
such as the fragment mass and charge distribution widths in heavy-ion collisions. In 
the limit where fluctuations around the TDHF path are small, then a good descrip- 
tion of these distributions is obtained within the time-dependent RPA (TDRPA) for- 
malism. The latter can be obtained from the Balian-Veneroni variational principle 
[109, 110]. Numerical applications to describe fragment mass and charge distribu- 
tions in deep-inelastic collisions have been recently performed and compared to ex- 
periment in Ref. [59] . Alternatively, such fluctuations could also be obtained from 
the stochastic mean-field (SMF) approach [111, 50], or from the time-dependent 
generator coordinate method (TDGCM) [112, 113]. 

Realistic calculations with these approaches beyond TDHF are often very de- 
manding from a numerical point of view and systematic applications are usually 
prohibitive even with modern high-performance computing facilities. In this chap- 
ter, we then focus on TDHF applications, paying attention to the limitations and 
possible improvements in the future. However, recent numerical calculations with 
the TDHFB and TDRPA approaches are also presented. 
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Nuclear molecules made of light nuclei such as C and O are formed and stud- 
ied thanks to nuclear collisions. The dynamics of the formation of such di-nuclear 
systems is the purpose of the present section. 



3.1 Structures infusion cross-sections 

Structures in fusion cross-sections are possible experimental signatures of nuclear 
molecules [1 14, 7]. However, structures in fusion excitation functions may also ap- 
pear in light systems which are not necessarily due to the formation of nuclear 
molecules. Such structures or oscillations appear clearly in cross-sections for the fu- 
sion of 12ch-I2c [115], 12ch-16o [116], and ^^O+^^O [117, 118, 119]. In particular, 
the discrete nature of angular momentum may reveal itself in fusion excitation func- 
tions as peaks associated to barriers for specific angular momenta [120, 121, 122]. 
The fusion cross-section is written as 

(£c.m. ) - L (2-^- + 1 )Pfus. [L, £c.m. ) (30) 

where pL is the reduced mass of the system. Pf^s.i^j^c.m.) is the fusion probability 
for the partial wave with orbital angular momentum L at the center-of-mass energy 

Ec.m. ■ 

TDHF calculations do not include tunnelling of the many-body wave-function, 
i.e., Pj,^^^^ = or 1. As a result, the fusion cross-section can be estimated with the 
quantum sharp cut-off formula [123] 

LmaxiEc.m) 

L=0 



2m£c., 



[W(£c■.m.) + l]^ (31) 



where L,nax{Ec.m.) is the maximum angular momentum at which fusion occurs at 
Ec.m.- For fusion of symmetric systems with 0+ ground-states, fusion can only occur 
for even values of the angular momentum. The cross-section with the sharp cut-off 
formula then reads 

0/H.V. {Ec.m. ) = [LmcLx {L,m,x + 3 ) + 2] . (32) 

An example of fusion cross-section in '^Oh-'^O as function of energy obtained 
with the TDHF3d code with the SLyAd Skyrme EDF [18] and Eq. (32) is shown 
in Fig. 1 with solid line. The sharp increases of the fusion cross-sections at the 
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Fig. 1 Fusion cross-section as 
a function of center-of-mass 
energy in 

16q^16q obtaingd 

with TDHF calculations. The 
cross-sections are computed 
with the sharp cut-off formula 
(solid line) and using Eq. (33) 
for the banier penetration 
probabilities (dashed line). 
The numbers indicate the 
position of the barriers B{L). 




25 
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positions of the angular momentum dependent barriers B{L) are due to the fact that 
fusion penetration probabilities are either or 1 at the TDHF level. 

These peaks are highly smoothen when tunnelling is taken into account. As a 
first approximation, one can estimate the barrier penetration probability according 
to the Hill-Wheeler formula [124] with a Fermi function 

P^„,.(L,£,,,,)~-^, (33) 

with XL = [E — B(L)]/e. Choosing the decay constant e = 0.4 [122], one gets the 
fusion cross-sections represented by a dashed line in Fig. 1 . Oscillations for Ec.m. > 
16 MeV are clearly visible and due to L-dependent barriers with L > \2Ti. Note that 
these oscillations are less visible for asymmetric systems due to the fact that all 
integer values of L are possible. In addition, the observation of these oscillations is 
limited to light systems up to, e.g., ^^Si-n^^Si [125, 122]. For heavier systems, the 
oscillations are indeed expected to be smeared out as the coupling to many reaction 
channels sets in [122]. 

To conclude, structures due to oscillations of the fusion cross-sections generated 
by the discrete nature of the angular momentum are expected to occur in light sys- 
tems, in particular for symmetric collisions. As a result, one should be careful in the 
search for resonances associated to molecular states in these systems. In particular, 
the observation of a peak in the fusion cross-sections may not be sufficient to as- 
sign such a structure to a resonance state. One should, in addition, search for this 
resonance in other channels, and investigate its decay properties [7]. 



3.2 Contact times around the barrier in ^^C+ °0 



Resonances at and below the barrier may strongly affect the reaction outcome. For 
instance, narrow resonances have been observed in the radiative capture '^C('^0,7)^**Si 
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close to the Coulomb barrier [7]. In particular, tentative spins of 4ti, 5ti and 6ti have 
been assigned in the collision energy range ~ 8.5 — 9 MeV [7, 1 14, 126]. 

To try to better understand in a dynamical way the presence of resonances and the 
contribution of relatively large spins at energies close to the Coulomb barrier, TDHF 
calculations have been performed on the '^C+'^O system [7]. For this system, the 
Coulomb barrier obtained with the TDHF3d code and the SLy4ii Skyrme functional 
[18] is B = 7.85 ±0.05 MeV. 

Fig. 2 shows the distance between the centers-of-mass of the fragments as a func- 
tion of time at Ec.m. =8.8 MeV for different values of the angular momentum. Fu- 
sion occurs at L <3h and re-separation at L > 4h. Without the presence of molec- 
ular states at 4h, 5h and 6h, the system would undergo a fast re-separation within 

1 — 2 zs. This indicates that these spins are populated by a direct transition toward 
resonant states of the compound nucleus. 

In the collision, the system may spend some time in a di-nuclear configuration, 
which presents some analogies with a molecular state. The excitation of the latter is 
then expected to increase with the lifetime of the di-nuclear system. The definition 
chosen here for the existence of a di-nucleus is the following: the nuclear density 
at the neck should be between 0.004 and 0.14 fm^^, i.e., lower than the saturation 
density of 0.16 fm^^. These nuclear densities correspond to distances between '^C 
and i''0 from 5.98 to 10.43 fm. 




The di-nuclear lifetime is shown as a function of the average angular momentum 
(L) in Fig. 3. Peaks are observed at the critical mean angular momentum {L)c for 
fusion. We get (L),. ~ 3.2^, 3.9h, and 4.4^1 for £,,,„. ^ 8.5, 8.8, and 9 MeV, respec- 
tively. Below (L)c , this time increases with (L) because the fusion process is slowed 
down by the centrifugal repulsion. In contrast, for (L) greater than (L)^, the frag- 
ments re-separate and the time of contact decreases for more peripheral collisions. 

With the present definition of the di-nuclear system lifetime, we see that such 
a system exists at relatively high angular momenta of 5h and bh during ^ 0.7 — 



Dynamical mean-field approach 



15 



xlO 



,-21 



1.8 



T 



Fig. 3 Di-nuclear lifetimes 
(see text) as a function of 
the mean angular momen- 
tum at E,.j„, = 8.5 MeV 
(blue circles), 8.8 MeV 




(green triangles), and 9 MeV 
(red squares). Adapted from 
Ref. [7]. 
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1.2 zs depending on the energy. The corresponding lifetimes are shorter than the 
Ufetimes of typical resonances at the Coulomb barrier, which are ^2.6 zs [127, 
126]. However, they may be sufficient to enable a direct excitation of a resonant 
or molecular state of the compound nucleus with a similar structure. In particular, 
we see that, for (L) = 5h and 6h (above the critical angular momenta), the lifetime 
increases with energy. Similarly, the probability for the population of a molecular 
state is also expected to increase. This may explain why the experimental cross- 
section for radiative capture is observed to increase with energy for these spins [7]. 



3.3 The J"" = 36+ resonance in ^^Mg+^^Mg 

The ^^MgH-^^Mg system presents a resonance with high spin {2>6Ti to 2>%h) at twice 
the Coulomb barrier, corresponding to an excitation energy E* ~ 60 MeV in the 
^^Cr At this energy, many decay channels are open and it is necessary, from an ex- 
perimental point of view, to investigate as many of these channels as possible. For 
instance, the = 36+ resonance in ^^MgH-^'*Mg is observed at Ec.m. = 45.7 MeV 
essentially in inelastic scattering [128] with a lifetime i = y — 3-9 zs. The de- 
cay of this resonance has been recently investigated in both inelastic and fusion- 
evaporation channels [8]. 

One particularity of the ^^Mg is its strong prolate deformation. Indeed, HF cal- 
culations with the EV8 code [129] and the SiLyAd Skyrme functional [18] of the 
Skyrme EDF [69] give a quadrupole deformation parameter ^2 — 0-4 [20]. As men- 
tioned in Sec. 2, one advantage of the TDHF theory is that it describes the struc- 
ture and the reaction mechanisms on the same footing. It is then well suited to 
investigate the role of the orientation of deformed nuclei on reaction mechanisms 
[20, 24, 29,31,52,53]. 
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In the present case, we can then investigate the role of the orientation of the 
^^Mg fragments at contact on the outcome of the reaction at L = 36h and fit .,,,. = 
45.7 MeV, that is, where the 36+ resonance is expected. In particular, we can deter- 
mine which relative orientations (e.g., tip-tip, tip-side, side-side) lead to the forma- 
tion of a di-nuclear system. 



Fig. 4 TDHF calculation of 
the reaction ^''Mg+^^'Mg 
at L = 36h and £[,„,. = 
45.7 MeV. The isoden- 
sity is plotted at half the 
saturation density, i.e., 
Po/2 = 0.08 fm-^. The defor- 
mation axes are aligned with 
the .t-axis at the initial time of 
the collision. 




Fig. 5 Same as Fig. 4 with 
the deformation axis of the 
right nucleus aligned with the 
y-axis at initial time. 




Fig. 6 Same as Fig. 4 with 
both deformation axes aligned 
with the v-axis at initial time. 




Figures 4, 5 and 6 present density evolutions of the ^'*MgH-^'*Mg system at L = 
36^ and Ec.,,,. = 45.7 MeV for different initial orientations of the nuclei computed 
with the TDHF3d code using the SLyAd Skyrme functional [18]. We observe that the 
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contact time between the fragments, if any, is extremely short as compared to the 
lifetime of the resonance, except when the contact occurs by the tip of the deformed 
collision partners, as it can be seen in Fig. 6. In the latter case, a di-nuclear system is 
formed with a lifetime similar or greater than the one of the resonance. As a result, 
it is likely to populate the resonance. This resonance could then be associated to 
a molecular state with the two fragments linked by a neck between their tips. In 
particular, it corresponds to a hyper-deformed nucleus as we can see in Fig. 6. 

To sum up, TDHF calculations provide an insight into the formation of the 36^ 
resonance in ^^Mg+^'*Mg. In particular, this resonance may be associated to a highly 
deformed state of the '*^Cr formed by the two aligned ^'*Mg in contact by their tips. 



4 a-clustering 

a-clusters play an important role in nuclear structure, both for ground- and excited 
states, due to the high binding energy of "^He [1, 2, 3, 4, 5, 130, 131, 132, 133]. They 
may also affect reactions. For instance, large cross-sections for the ^Li-n^'^^Pb^- 
2a+207Tl were recently measured [134], whereas this channel has a large negative 
Q-value. a-clustering in the entrance and/or exit channel is also believed to enhance 
a-transfer in heavy-ion collisions. One recent illustration is the formation of a- 
cluster states in ^'-Po with the reaction ^^^Pb('**0,''*C) [135]. a-clustering is also 
playing an important role in astrophysical processes such as helium burning in stars 
[136]. 

In this section, we present examples of reactions to illustrate the role of a-clusters 
in nuclear reactions. We first investigate the ^He-n^Be reaction of astrophysical in- 
terests, and, then, we study the survival of an a-cluster in its fusion with a ^''**Pb 
nucleus. 



4.1 Three- a cluster configurations in ^He+^Be 

The ^Be ground-state exhibits a 2-a-cluster configuration [1, 3]. This can be seen 
in Fig. 7 which shows the density profile of ^Be obtained from a HF calculation. In 
fact, the ^Be is unbound and decays by the emission of 2a with a lifetime of the 
order of ^ 10^'^ s. This lifetime, however, is long enough to allow the ^He-n^Be 
fusion by radiative capture to occur in stars. This reaction is considered as the main 
source of '^C production in the universe [136]. 

In this reaction, the '^C is formed in the state at E* = 7.654 MeV of excitation 
energy, i.e., just above the 3-a separation threshold. This state, also known as the 
Hoyle state [137], is believed to have a strong 3-a cluster configuration [1, 4, 5]. 
The exact configuration of the a-clusters in the Hoyle state is style under debate, in 
particular, it has been suggested that 3-a linear chains may contribute [138, 139]. 
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Fig. 7 HF ground-state den- 
sity of *Be obtained with the 
SLy4d parametrisation [18] of 
the Skyrme EDF. 




Fig. 8 Density snapshots of a 
*He-F'^Be reaction at £(■.„,, = 
2 MeV leading to a 3-0; linear 
chain (top) and evolving 
toward a triangular shape 
(middle) and then toward a 
more compact shape (bottom). 
Adapted from Ref. [38]. 
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Recently, Umar and collaborators have investigated the dynamics of 3-a linear 
chains formed in the ^He+^Be reaction using the TDHF formalism [38]. A 3-a 
linear chain is formed by the capture of a '*He by the tip of the ^Be in low-energy 
central collisions. Density profiles illustrating different stages of the dynamics of a 
3-a linear chain formed in such a central collision are shown in Fig. 8. The linear 
chain (top of Fig. 8) is present during a long time (up to 8 zs). During this time, 
the a-clusters present a longitudinal vibrational mode [38]. It is interesting to note 
that the clusters remain for such a long time, while the underlying formalism is a 
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mean-field model of independent particles, i.e., without imposing the presence of 
such clusters in the wave-function. 

For longer times, the chain becomes unstable due to the appearance of a bending 
motion favouring the formation of triangular shapes (middle of Fig. 8). The clusters 
still encounter some vibrations in this mode, with the center cluster oscillating per- 
pendicular to the left and right clusters. This vibration mode last for another 4 zs 
with almost no damping before a more contact shape is formed (bottom of Fig. 8). 

The role of impact parameter has also been studied by Umar et al. in Ref. [38]. A 
decrease of the lifetime of the linear chain with the impact parameter was observed, 
as reported in Fig. 9. We note that, although non-central collisions do not favour long 
lifetime of linear chains, they are stable enough to survive more than 1 zs as long 
as the collision occurs with the tip of ^Be and for impact parameters not exceeding 
0.5 fm. These times are of the same order of magnitude, if not larger, than typical 
lifetimes of di-nuclear systems formed in near-barrier heavy-ion collisions. 

To conclude, relatively long lifetimes of few zs are observed for linear chains of 
3-a clusters formed in "^HeH-^Be within the TDHF approach which does not assume 
a priori cluster components in the wave-function. The dynamics of these structures 
exhibit complex vibrational modes based on oscillations of the a-cores. In fact, 
similar vibrational modes have been found with the fermionic molecular dynamics 
approach by Furuta and collaborators [131]. In particular, possible strong effects on 
the vibrational response functions of light nuclei have been noticed. 



Fig. 9 Lifetime of the 3- 
a chain as a function of 
the impact parameter for 
the '^He+^Be reaction at 




0, 







0.1 



0.2 0.3 
Impact parameter (fm) 



0.4 



0.5 



= 2 MeV. Adapted from 
Ref [38]. 
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4.2 Survival of a-clusters in ^He+^'^^Pb near-barrier fusion 



The previous section emphasises the survival of a-clustering after a capture process 
in a light system, a-clustering is not limited, however, to light nuclei. Indeed, the 
well-known a-radioactivity, which occurs essentially in heavy nuclei, is another 
form of a-clustering. In addition, the reverse process of capture of an a by a heavy 
nucleus in a transfer or a fusion reaction may form excited states interpreted as 
nucleusH-a molecules. For instance, new excited states have been recently observed 
in ^'^Po that are interpreted as a+^^^^Pb configurations [135]. 

The dynamics of the a+^*^*Pb system after capture of the a by the heavy partner 
has been investigated with an early TDHF code in Ref. [14]. However, these cal- 
culations were performed with a simplified Skyrme functional. In particular, they 
did not include the spin-orbit interaction which is known to be crucial for a proper 
description of nuclear reactions [140, 21, 22]. In fact, the spin-orbit interaction and 
the difference between proton and neutron mean-fields may induce a "dissolution" 
of an a-particle entering the mean-field of a collision partner [141]. 

To test this idea, three-dimensional TDHF calculations with the TDHF3d code 
and the SLy4c/ Skyrme functional [18] including spin-orbit terms have been per- 
formed. The ^He-H^^^Pb central collision at Ec.,,,. = 19.8 MeV, i.e., slightly above 
the Coulomb barrier, is considered. The evolution of the density associated to the 
single-particle wave-functions initially forming the "^He projectile is shown with 
green contours in Fig. 10. A clear dissolution of the a-particle in at least two compo- 
nents is observed as soon as it enters the mean-field of the heavy partner, represented 
by a single isodensity at half the saturation density (white solid line in Fig. 10). 

To conclude, the correlations responsible for the survival of a-clustering in the 
previous study of the ^He-n^Be system are not strong enough, at the mean-field level, 
to enable similar effects in a ^He-n^^^Pb central collision at the barrier. Indeed, the 
2p-2n spatial correlations in the ''^He fragment are lost in the mean-field of the heavy 
partner. The experimentally observed a-cluster states in the ^'^Po nucleus [135] 
may either be due to beyond-mean-field correlations, or to an a-cluster located at 
the surface and orbiting around the heavy core. To test the latter idea, the present 
calculations should be repeated above the barrier and around the grazing angle to 
investigate possible long-lived a-cluster configurations. 



5 Transfer in heavy-ion collisions 

Transfer reactions are highly sensitive to cluster effects. This is illustrated in the 
a-transfer experiment of Ref. [135] discussed above. Other types of clusters may 
also be transferred in heavy-ion collisions, such as pairs of protons and neutrons 
[142]. These types of multi-nucleon transfer reactions are in competition with the 
sequential transfer, i.e., a transfer of independent nucleons. 

Theoretically, transfer probabilities in microscopic approaches are obtained from 
a particle number projection technique. First, this technique is described. It is 





Fig. 10 Snapshots of the density profile in the "^He+^'^^Pb central collision at E^.,,,,. = 19.8 MeV. 
Wave-functions initially belonging to the '^He are shown with green contours. A single isodensity 
at half the saturation density (white solid line) is shown for the wave-functions initially belonging 
to the ^"*Pb fragment. 
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then applied at the TDHF level to estimate sequential transfer probabilities in the 
i6Q^208p|^ system. The latter are used as a reference to interpret experimental data 
in terms of cluster transfer. Finally, pairing vibrations, which may be excited in pair 
transfer reactions, are studied with the TDHFB approach. 



5.1 Particle number projection technique 

Transfer probabilities have been extracted at the TDHF level [143, 56] thanks to 
a projection onto a good particle number technique. This technique is standard in 
beyond-mean-field models for nuclear structure subject to pairing correlations [87, 
144]. Here, it is applied on the outgoing fragments to determine their proton and 
neutron number probabilities. 

It is possible to extract the component of the wave function associated to a spe- 
cific transfer channel using a particle number projector onto protons or neutrons 
in the x > region where one fragment is located at the final time, the other one 
being in the .y < region. Such a projector is written [56] 

pRiN) ^ ^ TdO (34) 

where 

Nr=Y, [dr a\rs)a{rs)0{x) (35) 

counts the number of particles in the x > Q region (0 (x) = 1 if x > and (x) = 
elsewhere). Isospin is omitted to simplify the notation. 

The projector defined in Eq. (34) can be used to compute the probability to find 

nucleons in x > in the final state |^), 

9 1 r^'' 

\PR{Nmr^:^ de e-'''{^Me)), (36) 
2n Jq 

where \^r{9)) = e'^^^lif)) represents a rotation of |0) by a gauge angle 9 in the 
gauge space associated to the particle number degree of freedom in x > 0. Note 
that is an independent particle state. The last term in Eq. (36) is then the 

determinant of the matrix of the occupied single particle state overlaps [56]: 

{(^\MO))=detiF) (37) 

with 

Fij=Y, /dr(pr(r)(p;(ry^®«. (38) 

The integral in Eq. (36) is discretised using 9„ = 2 TTn/M with the integern = 1 ■■■M. 
Choosing M ~ 300 ensures numerical convergence for the '^OH-^''^Pb system. 
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We mentioned in Sec. 4 the a-transfer reaction ^*'^Pb('^0,'^C) populating a-cluster 
states in ^'^Po [135]. One could expect a similar a-transfer mechanism to dominate 
the 208pb(i6o,'*C) reaction channels below and around the barrier [145, 146, 147]. 
Indeed, '^^O is a good candidate for a-clustering, and a-condensates [148] as well as 
linear a-chains [132] have been predicted. Experimental indications of a-clustering 
in '^O have also been reported [149, 150]. Recent experimental data [57] showed, 
however, that the most probable carbon isotope formed in *^OH-^*^^Pb is '''^C, indi- 
cating a dominance of two-proton transfer against a-transfer 



Fig. 11 Charge transfer prob- 
ability in I'-O+^o^Pb. TDHF 
calculations of the AZ — 1 
(Ip stripping) and AZ = 2 
(2p stripping) transfer prob- 
abilities as a function of the 
distance of closest approach 
r„,i„ are shown by the full 
curves. Experimental data are 
shown for AZ = 1 (squares) 
and AZ = 2 (diamonds). The 
vertical dashed line indicates 
the average barrier radius. 
The dotted curve shows the 
square of the TDHF AZ = 1 
transfer probabilities. The red 
dashed curve shows the sub- 
barrier fusion probability as 
calculated with the coupled- 
channels code CCFULL. 
From Ref. [57]. 
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TDHF calculations [56] have been performed with the TDHF3d code and the 
Shy Ad Skyrme functional [ 1 8] to estimate the probability for production of nitrogen 
{AZ = 1) and carbon {AZ = 2) in ^^0+^^^¥h. The results are shown in Fig. 1 1 (sohd 
lines) as a function of the distance of closest approach for Coulomb trajectories 

We observe that these probabilities decrease exponentially with increasing r,„,„, as 
expected for quasi-elastic transfer [151]. 

Due to the independent particle approximation, the TDHF approach is only able 
to describe sequential multi-nucleon transfer, i.e., neglecting cluster correlations. 
The probability for sequential transfer of two nucleons is sometimes approximated 
by the square of the one-nucleon transfer probability. This leads, however, to an 
overestimation of the two-proton sequential transfer in ^^O-n^'^'^Pb, as can be seen in 



-cosec^^ I . (39) 
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Fig. 1 1 (compare dotted and solid black lines). In fact, the above approximation is 
valid when a large number of particles are available for transfer toward states with 
large degeneracies. In general, this criterion is not fulfilled due essentially to the 
relatively small number of single-particle states around the Fermi level. Indeed, the 
latter have the smallest binding energies, and, then, the largest transfer probabili- 
ties. In the '^O case, for instance, transfer is dominated by the lpi/2 single particle 
states. Quantum microscopic approaches such as the TDHF theory are then needed 
to estimate correct sequential multi-nucleon transfer probabilities. 

Comparison with experimental data in Fig. 1 1 indicates that to assume pure se- 
quential transfer, as in TDHF calculations, leads to an overestimation of one-proton 
transfer probabilities by a factor ^ 2, and to an underestimation of two-proton trans- 
fer probabilities by approximatively one-order of magnitude. This is interpreted as a 
strong cluster effect in "^O-n^^'^Pb charge-transfer [56, 57]. Indeed, correlations such 
as proton pairing or a-clustering favour the transfer of two protons, while they re- 
duce the probability for transferring only one proton. In fact, it is shown in Ref. [57] 
that pairing correlations dominate the AZ = 2 channel over a-clustering. 

It is interesting to note that at large distances of closest approach, i.e., r,„„, > 
13 fm, the sum of the AZ ~ I and 2 channels is rather well reproduced by the 
TDHF calculations [56]. Closer to the barrier, however, sub-barrier fusion domi- 
nates, as shown by the coupled-channels calculations with the CCFULL code [152]. 
As a result, the experimental transfer probabilities are reduced at the barrier, induc- 
ing a deviation from the exponential dependence observed at larger distances. One 
drawback of the TDHF approach is that it does not enable tunnelling of the many- 
body wave-function. Thus, sub-barrier fusion is not included in TDHF calculations, 
inducing an overestimation of the total transfer probabilities close to (but below) the 
barrier 

To sum up, the TDHF approach provides a good estimate of sequential transfer 
probabilities well below the barrier. These sequential transfer probabilities can be 
used as a benchmark to compare with experimental data in order to emphasise the 
role of cluster correlations on the transfer mechanism. 



5.3 Pairing vibrations 

It is possible to include pairing correlations in the mean-field dynamics by con- 
sidering quasiparticle vacua instead of independent particle states. This leads to 
the TD-BCS model when only pairs between time-reversed single-particle states 
are considered, or, more generally, to the time-dependent Hartree-Fock-Bogoliubov 
(TDHFB) formalism [153]. Numerical applications are now possible thanks to the 
recent development of TD-BCS [88, 89] and TDHFB [90,91, 80] codes. 

A proper description of pairing dynamics is crucial to investigate the evolution 
of nuclei produced by a pair transfer mechanism [154, 142, 155, 156, 157, 158]. 
Indeed, pair transfer reactions are a good tool to excite the so-called "pairing vi- 
brations" [159, 87, 160, 154, 142]. Pairing correlations are then expected to induce 
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a collectivity which manifests itself as an increase of transition amplitude toward 
these states. 

Recent studies of pairing vibrations have been performed at the QRPA level 
[155, 156], that is, in the linearised version of TDHFB [153, 87]. Here, we dis- 
cuss a similar study with a real-time description of pairing vibrations excited in 
two-neutron transfer reactions in ^^Ca with a fully self-consistent TDHFB code 
[90, 161, 162]. Applications to other systems including oxygen, calcium, and tin 
isotopes can be found in Refs. [90, 162, 163]. 

Starting with an even-even nucleus ground-state with A nucleons and spin-parity 
0+, and assuming a AL = direct pair transfer reaction, pair vibration states with 

= 0+ are populated in the A + 2 (pair addition) and/or A — 2 (pair removal) nuclei. 
The transfer process is simulated by an initial excitation generated by a boost with 
a Hermitean pair-transfer operator [160] 



where the arrows label the spin of the single-particles (we omit the isospin to sim- 
plify the notation). In the present application, f{r) is a Fermi-Dirac spatial distri- 
bution containing the nucleus and cutting at 4 fm outside the nucleus. Its role is to 
remove unphysical high energy modes associated to pair creation outside the nu- 
cleus. 

The evolution of the system after the boost is obtained with the TDHFB RAD 
code [90]. This code solves the TDHFB equation in spherical symmetry with a 
full Skyrme EDF and density-dependent pairing effective interaction. The linear 
response of {F){t) after a boost excitation is shown in Fig. 12 for a '^'^Ca ini- 
tial ground-state. The SLy4 parametrisation of the Skyrme EDF is used [74], to- 
gether with a "volume" pairing effective interaction of the form fo5(f i — fo) with 
fo = — 187 MeV.fm^^ and a cut-off energy of 80 MeV in the quasi-particle spectrum 
to avoid ultra-violet divergence. See Refs. [163, 164, 162] for a discussion on the 
role of the form of the pairing functional on the excitation of pairing vibrations. 

The apparent chaotic behaviour of (F) (r) in Fig. 12 is due to the simultaneous ex- 
citation of several states. A simple Fourier analysis can be used to extract the energy 
and relative contributions of these states to the time evolution of {F){t). The result- 
ing strength function is plotted in Fig. 13 with a solid line. Both pair additional and 
pair removal (indicated by the arrows) modes are present. The unperturbed strength 
function (dashed line) obtained by removing the self-consistency of the generalised 
mean-field is also shown. Overall, an increase of the strength is observed due to 
the dynamical pairing residual interaction present in TDHFB, but neglected in the 
unperturbed response. This increase of the strength is a signature for collective mo- 
tion, indicating that several quasi-particles participate to the vibrational modes. In 
addition, this residual interaction lowers the transition energies due to its attractive 
nature. These characteristics of pairing vibrations are in agreement with previous 
observations with the QRPA model [155]. 

To conclude, collective pairing vibrations may be excited in pair-transfer reac- 
tions. These reactions are then a unique probe to investigate the dynamics of pairs 
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Fig. 12 Evolution of {F){t) 
after a pair transfer type 
excitation on '^''Ca. The inset 
shows the same quantity at 
early times. Adapted from 
Ref. [163]. 
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Fig. 13 TDHFB strength 
function (solid line) associ- 
ated to the evolution of (F)(f) 
in Fig. 12. The unperturbed 
spectnam (dotted line) is also 
shown. The arrows indicate 
pair removal transitions from 
deep hole states. Adapted 
from Ref [163]. 
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of nucleons in nuclei, in particular their vibrational modes. Pairing vibrations have 
been studied at the TDHFB level with spherical symmetry. The recent development 
of a three-dimensional TDHFB code [91] might enable similar studies in deformed 
nuclei and for L / modes. The complete description of pairing vibration excita- 
tions in heavy-ion collisions might also be possible with such a code. 



6 Deep-inelastic collisions 

Deep-inelastic collisions (DIC) have been widely studied in the past [165, 166, 167, 
168, 169]. They are characterised by a strong damping of the relative kinetic energy 
and an orbiting of the di-nuclear system before re-separation of the two fragments. 
In particular, large widths of the fragment mass and charge distributions are usually 
observed. 
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Die have been used to investigate isospin equilibration in damped collisions 
[170, 171, 172, 173], and to produce nuclei and study their structure (see for ex- 
ample Refs. [174, 175]). Upcoming radioactive beams will be used to further inves- 
tigate transport properties of isospin asymmetric nuclear matter [176]. 

The characteristics of DIC provide stringent tests to nuclear transport models 
[177, 178, 179]. For instance, it has been shown that standard TDHF calculations 
usually fail to reproduce the large widths of mass and charge distributions [143, 
177]. This is due to the fact that the many -body wave-function is constrained to be 
a single Slater determinant at all time [180]. Fluctuations are then computed with 
beyond TDHF approaches such as the TDRPA [110] and stochastic mean-field [111] 
formalisms. 

We first discuss briefly the calculation of fluctuations at the TDRPA level. Then, 
we present applications to the '*''CaH-'*''Ca reaction well above the barrier. 



6.1 Fluctuations of one-body observables 

Balian and Veneroni (BV) have introduced a variational principle in which the 
TDHF theory turns out to be optimised to the expectation value of one-body ob- 
servables [109]. It could then fail to reproduce quantities like two-body observables 
and fluctuations of one-body observables. 

Balian and Veneroni also used their variational principle to derive a prescription 
for fluctuations and correlations between one-body observables [110, 181] (a de- 
tailed derivation can also be found in Ref. [13]). This prescription is, in fact, fully 
equivalent to the TDRPA approach where small fluctuations around the mean-field 
evolution are considered. 

The BV variational principle can then be used to determine an optimum pre- 
diction for correlations Oxy and fluctuations Oxx of one-body operators assuming 
small fluctuations around the mean-field path [110, 181]. Correlations are obtained 
from the general expression 



where Gxy has the sign of the term between the absolute value bars. The X and 
Y distributions are correlated (resp. anti-correlated) for Oxy > (resp. Gxy < 0). 
Fluctuations are obtained by taking X ^Y, leading to 



Assuming independent particle states, the BV variational principle leads to the 
prescription 




-{{XY) + {YX})-{X){Y) 



(41) 




(42) 



a|y (fi ) = Hm — 2 tr{ [p (to) - px (fo, e)] [p (to) - p}'(fo, e)] }, 



(43) 
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where tr denotes a trace in the single-particle space. The one-body density matrices 
px (f , e ) obey the TDHF equation ( 1 ) with the boundary condition at the final time f i 

Px(fl,e) = e'^''^p(rl)e-'■^'?^ (44) 

while p(f) = px(f ,0) is the standard TDHF solution. 

Eq. (43) has been solved numerically in the past for particle number fluctua- 
tions with simple effective interactions and geometry restrictions [182, 183, 184]. 
Modern three-dimensional TDHF codes with full Skyrme functionals are now used 
for realistic applications of the BV variational principle to determine these fluctu- 
ations [185, 186, 59, 187, 13] as well as the correlations between the proton and 
neutron number distributions [59, 187]. See Ref. [13] for numerical details of the 
implementation of Eq. (43). 



6.2 The ^^Ca+^'^Ca reaction well above the barrier 




The Can- Ca reaction has been investigated at Ec.m. = 128 MeV (^^ 2.4 times 
the barrier) [59] with the TDHF3d code and its TDRPA extension using the SLy4^f 
Skyrme functional [18]. Fig. 14 provides some examples of trajectories obtained 
at different angular momenta. The TDHF calculations for this reaction predict that 
fusion occurs at L < 66 [59]. We see in Fig. 14 that orbiting followed by re- 
separation is predicted at L = 7Qh. Partial orbiting at L = 15h and smaller nuclear 
deflections at larger L are also observed. 

It is interesting to note that different angular momenta may lead to similar scat- 
tering angles. This is the case, for instance, with L = lOTi and L = 15Ti which are 
associated to different orbiting trajectories (see Fig. 14). In fact, the amount of or- 
biting changes very rapidly with L for DIC. Comparisons with experimental data 
imply then to perform calculations with a small angular momentum step in the or- 
biting region [59]. 

As mentioned before, DIC are not only characterised by a large orbiting of the 
di-nuclear system. They are also associated to a large damping of the initial relative 
kinetic energy. This is quantified by the total kinetic energy loss TKEL ~ Ecm. — 
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El — E2, where Ei 2 are the asymptotic kinetic energies of the fragments in the exit 
channel. The TKEL in '*°Ca+'*"Ca at £,,„,. = 128 MeV ai-e shown in Fig. 15 as a 
function of the initial angular momentum. The maximum TKEL of ~ 60 — 70 MeV 
are obtained close to the critical angular momentum for fusion. These TKEL have to 
be compared with the Viola systematics for fission fragments [188]. The latter gives 
an expected TKEL in symmetric fission of ^76 MeV. This indicates that the DIC 
around L ~ IQTi are almost fully damped. Note that this result is obtained with TDHF 
calculations which contain one-body dissipation only. As a result, the damping of 
relative kinetic energy in DIC is essentially of one-body nature. 

Another characteristic of DIC is the large width of the fragment mass and charge 
distributions. Independent particle descriptions such as the TDHF theory usually 
strongly underestimate these widths [143, 177, 59]. It is then necessary to include 
beyond TDHF fluctuations with the TDRPA [110] or with the stochastic mean-field 
(SMF) approach [111]. Calculations with the TDRPA [182, 183, 184, 185, 186, 59, 
187] and with the SMF approach [50, 58] indeed predict larger fluctuations than the 
TDHF flieory. 

Fig. 16 Fragment charge 
fluctuation Gzz with TDHF 
(dashed line) and TDRPA 
(filled circles) as a function 
of the angular momentum L 
in units of h. The fluctuations 
Gnn of the neutron number 
distributions, not shown, are 
very close to the proton ones. 
Correlations Gmz between 
proton and neutron num- 
bers distributions obtained 
with TDRPA calculations are 
shown with open triangles. 
Adapted from Ref. [59]. 
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This is illustrated in Fig. 16 where the fragment charge fluctuations obtained 
from TDHF (dashed line) and TDRPA (filled circles) are reported as a function 
of the angular momentum for the '*''Ca+'**'Ca reaction at Ecm. = 128 MeV [59]. 
We see that, for large angular momenta L > 9QTi, both approaches predict similar 
fluctuations. These reactions are very peripheral and associated to small TKEL of 
few MeV (see Fig. 15). Collisions at L > 90/i are then dominated by quasi-elastic 
scattering. This shows that the TDHF approach may be used safely to compute 
quasi-elastic transfer (see also section 5) as it provides similar fluctuations than the 
TDRPA for these quasi-elastic events. 



Fig. 17 Fluctuations Ozz of 
the proton number distribu- 
tions obtained with TDHF 
(dashed line) and TDRPA 
(solid line) calculations are 
plotted as function of the 
center-of-mass scattering an- 
gle 9c.,„.- Experimental data 
(open circles) from Ref. [165] 
are also indicated. Adapted 
from Ref [59]. 
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On the other hand, we observe in Fig. 16 that fluctuations with the TDRPA are 
much larger than with the TDHF approach for more central collisions, in particular 
in the DIC region at L '-^ 7Qh. A comparison of these fluctuations with the experi- 
mental data of Roynette and collaborators [165] has been performed in Ref. [59]. 
The results are reported in Fig. 17 for a selection of events with TKEL > 30 IVleV. 
Although the TDRPA results still underestimate experimental data, they provide 
both a better qualitative and quantitative agreement than the TDHF calculations. In 
fact, the plateau observed at large angles contains a contribution from fusion-fission 
events [165]. The latter are not treated in the calculations and may be the origin 
of the remaining difference between the TDRPA calculations and the experimental 
data [59]. 

Correlations between proton and neutron numbers distributions have also been 
computed recently with the TDRPA approach for '^^Ca-n^^'Ca collisions at E^j,,, = 
128 MeV [59]. In standard TDHF calculations, these correlations are strictly zero. 
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In the TDRPA, however, they become important in the DIC region, as shown in 
Fig. 16 (open triangles). Although they are negligible for quasi-elastic scattering 
(L > 9Qh), they are similar to the charge fluctuations for the most damped events. 
This indicates that, e.g., an addition of several protons in one fragment is likely to 
be accompanied by an addition of neutrons as well. This is a manifestation of the 
symmetry energy which favours N = Z fragments. 

Finally, combining neutron and proton fluctuations with their correlations, one 
can estimate the distribution of nuclei produced in the reaction. Let us assume a 
Gaussian probability distribution of the form 



^(n,z) = ^(0,0)exp 



1 



« z 2pnz \ 



(45) 



where n and z are the number of transferred neutrons and protons, respectively. The 
probability for the inelastic channels reads 



),0)= ( 27t(7N(7z\/ I - p- 

The dimensionless quantity 



:sign(cT^z)-^= /"'^ , (46) 



quantifies the correlations and obeys |p| < 1. The case p —0 corresponds to inde- 
pendent distributions of the form £P{n,z) = 3^{n)3^{z). On the other side, fully 
(anti-)correlated distributions are found in the limit p — > ± 1 . 

Iso-probabilities corresponding to 3^{n,z) = ^{0,0)/2 are plotted in Fig. 18 
for different angular momenta. Independent proton and neutron distributions would 
produce ellipses with the principal axis parallel to the abscissa or to the ordinate. 
This is not what is observed, particularly for the smallest angular momenta corre- 
sponding to the most violent collisions. We see that not only the fluctuations are 
important to determine distributions of DIC, but the correlations play a significant 
role as well. 

To conclude, the TDHF theory is a good tool to compute transfer probabilities 
in quasi-elastic scattering. However, beyond TDHF fluctuations are mandatory to 
describe fragment mass and charge distributions in more violent reactions such as 
deep-inelastic collisions. Calculations based on the TDRPA indeed provide a better 
agreement with experimental data than standard TDHF codes. The correlations are 
also shown to be important in DIC with these TDRPA calculations. They should be 
sensitive to the symmetry energy and might be used in the future to test this part of 
the functional, in particular with exotic beams. 
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7 The quasi-fission process 

When two nuclei collide with an energy above the Coulomb barrier and a small 
enough impact parameter, a capture of the two fragments is expected to occur, i.e., 
a di-nuclear system is formed after dissipation of the relative kinetic energy [9]. 
The outcome of such a di-nuclear system is either (/) fusion, i.e., the formation 
of a unique system where two-centers cannot be identified anymore in the density 
distribution, or (//) a re-separation after a possible multi-nucleon transfer between 
the fragments. 

In light and medium mass systems, fusion is usually enabled by a close contact 
between the fragments. For heavy systems with typical charge products greater than 
1600- 1800, however, the second process is often dominant around the Coulomb 
barrier, leading to a fusion hindrance in these systems [189]. Instead of fusing, 
the di-nuclear system encounters a re-separation in two fragments after a possi- 
ble exchange of a large number of nucleons. This process is called quasi-fission 
as the characteristics of the fragments may exhibit some strong similarities with 
those emitted in statistical fission of the compound nucleus formed by fusion 
[190, 191, 192]. Note that, although much less probable than fusion, quasi-fission 
may also occur in lighter systems [193, 194, 195]. 

Firstly, we present some TDHF calculations of fusion hindrance in several heavy 
systems. Then, we investigate the effect of some structure properties of the collision 
partners, in particular their deformation, on the quasi-fission process. 
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7.1 Fusion hindrance in heavy systems 

Let us illustrate the phenomenon of fusion hindrance with TDHF calculations of 
heavy systems using the TDHF3d code with the SLy4c/ Skyrme functional [18]. 
Fig. 19 shows the density evolution of a ^^Fe+^^^Pb (Z1Z2 = 2132) central collision 
at Ec.in. = 240 MeV. This energy is well above the Coulomb barrier Indeed, the 
barrier computed with the proximity model [ 1 96] is — 224 MeV. However, this 
collision does not lead to fusion. Indeed, despite the formation of a rather compact 
di-nuclear system, the latter end up in quasi-fission. The lifetime of the di-nuclear 
system is 3 zs, which is too short to enable a full mass equilibration [190, 191, 
192, 197] which would be signed by a symmetric exit channel. 



Fig. 19 Snapshots of the den- 
sity for a 56Fe-l-2"*'Pb central 
collision at Ec.,,,. = 240 MeV. 
The isodensity is plotted at 
half the saturation density, 
i.e., at po/2 = 0.08 fm^^. 
Snapshots are shown every 
0.75 zs. Time runs from top to 
bottom. 




Fusion may eventually occur in some collisions if enough energy above the 
Coulomb barrier is brought into the system. This additional energy, sometimes 
called "extra-push" energy, may be computed with phenomenological approaches 
such as the extra-push model of Swiatecki [198]. Modern TDHF calculations are 
also able to determine such fusion thresholds without any parameter adjusted on 
reaction mechanism [13, 47, 199]. 

An example is provided in Fig. 20 where the neck density is plotted as a function 
of time for ^"Zr-n'^^Sn (Z1Z2 = 2000) central collisions at several energies. This 
system has a proximity barrier Bprox. — 215 MeV. Densities exceeding 0.13 fm^^ 
are observed at this energy and above. However, re-separation occurs for energies 
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Fig. 20 Neck density in 
^Zr+'^'^Sn central collisions 
as function of time [200]. The 
center-of-mass energies are 
indicated in the legend. 
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smaller than 240 MeV. Increasing contact times with energy are observed below 
240 MeV. At Ee.m. = 240 MeV, however, the neck survives more than 14 zs, which 
may be interpreted as a fusion process. An extra-push energy of 22 to 25 MeV above 
the proximity barrier is then needed for the ^''Zr-i-'^^Sn system to fuse. 



Fig. 21 Distance between the 
centers-of-mass of the frag- 
ment in ^"Zn+-"'*Pb central 
collisions as a function of 
time for center-of-mass ener- 
gies between 240 MeV and 
380 MeV [200]. 
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Increasing the collision energy does not guarantee to reach such a fusion thresh- 
old in all systems. For instance, TDHF calculations predict that fusion is not possible 
in the ™ZnH-2''^Pb (Z1Z2 = 2460) system [48]. This is illustrated in Fig. 21 where 
the distance between the centers-of-mass of the fragments is plotted as a function of 
time at different center-of-mass energies ranging from 240 MeV to 380 MeV. The 
proximity barrier for this system is Bprox. — 252 MeV. We see that, at an energy of 
more than 100 MeV above this barrier, the system is still not fusing. In fact, we 
observe a rise and fall of the contact time in this system with increasing energy in 
Fig. 22. This indicates that no fusion is expected in this system. 
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Fig. 22 Contact time between 
the fragments arbitrarily de- 
fined as the time the systems 
spend at a distance between 
their centers-of-mass smaller 
than 15 fm. Adapted from 




.3 



Ref. [48]. The arrow indicates g 
a lower limit. c.m. prox. 



To conclude, the fusion hindrance in heavy systems which hinders the formation 
of very heavy elements by fusion is described in the TDHF approach. When not 
fusing, the di-nuclear systems encounter quasi-fission within several zeptoseconds. 
Fusion may occur in some systems with an additional extra-push energy. However, 
others, with larger charge products are never found to fuse whatever the energy. 



7.2 Effects of the structure of the collision partners 

The previous section shows the importance of quasi-fission in the outcome of heavy 
di-nuclear systems. Realistic descriptions of quasi-fission are challenging because 
many degrees of freedom are at play. In addition, the shape of the di-nuclear system 
and its mass and isospin repartition evolves dynamically on different time scales. 

The equilibration of the isospin degree of freedom has been studied with the 
TDHF approach in several works [201, 19, 27, 28, 51, 54, 55, 60, 44, 59, 61]. It has 
been shown to occur on a typical time scale smaller than 2 zs [61]. This time scale is 
smaller than standard quasi-fission times. Quasi-fission fragments are then expected 
to have a similar A^/Z ratio than the compound nucleus. 

Shell effects of the collision partners have also been shown to play an important 
role in the competition between fusion and quasi-fission [202, 203, 61]. In particular, 
fusion might be eased by the magicity as less dissipation is expected in the fusion 
valley of magic nuclei [204, 205]. As a result, more compact di-nuclear systems 
can be formed [206, 207, 203, 208]. This may explain the success of super-heavy 
element synthesis with the doubly-magic ^^Ca projectile [209, 210]. However, these 
effects remain to be investigated with quantum microscopic approaches such as the 
TDHF formalism. 

The role of deformation and orientation on quasi-fission has been investigated in 
several experiments [211, 212, 193, 213, 214, 215, 216]. These studies led to the 
general conclusion that quasi-fission in collisions with the tip of a deformed heavy 
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target (e.g., in the actinide region) is dominant. On the other side, collisions with 
the side lead to more compact shapes which favour long lifetimes of the di-nuclear 
systems and then increase the fusion probability. 

The effect of deformation and orientation on quasi-fission has been investigated 
recently with TDHF calculations [13, 47, 48, 217] of the '^''Ca+^^^U system. In 
particular, it has been shown that collisions with the tip of ^^^U do not lead to fusion, 
but to quasi-fission with a contact time of 10 zs and a partial mass equilibration 
almost independent with energy. On the contrary, contact with the side of ^^^U leads 
to long di-nuclear lifetimes above the barrier which may induce the formation of a 
compound nucleus. 

In Ref. [61], it was shown that ^"Ca and '^'^Ca behave differently as far as the 
interplay between quasi-fission and fusion is concerned when colliding with a ^"^^Pb 
target. In particular, the hindrance of quasi-fission due to shell effects is observed 
only with ^^Ca, despite the fact that ^''Ca is also doubly magic. This is interpreted 
as an effect of isospin asymmetry which, in the case of ^''Ca, induces a fast N jZ 
equilibration breaking the magicity of the fragments in the di-nuclear system [61]. 



Fig. 23 Contact time be- 
tween the fragments as 
function of center-of-mass 
energy normalised by the 
proximity barrier [196] in 

40.48 ca+238u 

sions. Here, contact times 
are defined as the time dur- 
ing which the neck density 
exceeds half the saturation 
density po/2 = 0.08 fm^l 
Collisions with the tip (solid 
line) and the side (dashed line) 
of the ^^'^U are considered. 
The arrows indicate lower 
limits. 
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It is interesting to see if differences in the quasi-fission process between the two 
calcium isotopes are also observed in collision with a heavy deformed nucleus. 
Contact times, defined, here, as the time during which the neck density exceeds 
Po/2 = 0.08 fm^^, have been computed in '*^CaH-^^^U central collisions with the 
TDHF3d code and the SLy4fif Skyrme functional [18]. The evolution of these con- 
tact times as a function of energy is plotted in Fig. 23 for collisions with the tip 
(solid line) and with the side (dashed line) of ^^**U. The behaviours of the contact 
times present similarities between the two isotopes. For instance, quasi-fission times 
for collisions with the tip are of the order of 10 zs for both '^'^Ca and ^^Ca. Collisions 
with the side also present a sharp increase of the contact time above the barrier in 
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both cases, with long contact times (more than 20 zs) at high energy which could 
lead to fusion. However, long contact times for collisions with the side are reached 
at higher energies with "^^Ca than with ^"Ca. This might be attributed to shell effects 
in the ^^Ca-like fragment which are absent in the '^''Ca-like fragment due X.o N /Z 
equilibration [61]. This effect needs further investigations. 

To sum up, the competition between quasi-fission and fusion is affected by the 
structure of the collision partners. In particular, the deformation and the orientation 
is crucial. Recent TDHF calculations confirm that collisions with the tip lead essen- 
tially to quasi-fission, while long contact times possibly leading to fusion may be 
reached above the barrier for collisions with the side producing to more compact 
shapes. More investigations with quantum microscopic approaches are needed to 
gain a deeper understanding on how the various structure characteristics, such as, 
e.g., shell effects and isospin, affect quasi-fission. 



8 Actinide collisions 

Collisions of actinides form "hyper-heavy molecules" with ~ 500 nucleons in inter- 
action during short times of few zeptoseconds [10]. The description of their dynam- 
ics is of course a great challenge for theorists. These reactions may be an alternative 
way to produce more neutron-rich heavy and super-heavy nuclei than those formed 
by fusion [218, 219, 53, 187]. This is possible thanks to the fact that actinides have 
large neutron to proton ratio, of the order ofN/Z^ 1 .5. 

Another interest of actinide collisions is the possibility to make the QED vacuum 
unstable due to the strong electric field [220, 22 1 , 222] . As a result, a spontaneous 
decay of the vacuum by the emission of a e+e^ pair is expected. The lifetime of the 
hyper-heavy molecule is a crucial parameter which determines if this QED vacuum 
decay may be observed experimentally. 

These applications of actinide collisions require a precise description of the 
reaction mechanisms. The dynamics of actinide collisions has been investigated 
with various theoretical approaches, including macroscopic models [219, 223, 224], 
semi-classical microscopic approaches [225, 226], and the TDHF theory [16, 17, 
52,53]. 

In the following we describe, first, the role of the relative orientation of the nu- 
clei on the reaction mechanisms. Then, we look for the optimal conditions for the 
observation of spontaneous e+e^ emission. 



8.1 Role of the initial orientation 

Di-nuclear systems formed in actinide collisions are too heavy to fuse. They always 
encounter quasi-fission. As in quasi-fission with lighter projectiles, the deforma- 
tion and orientation of the nuclei play a crucial role in the outcome of the collision 
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[10]. This has been confirmed with recent TDHF calculations of the 238jj^238jj ^^2] 
and ^^^Th+^^^Cf [53] systems with the TDHF3d code and the SLy4(i Skyrme func- 
tional [18]. 
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Fig. 24 Isodensities at half the saturation density, i.e., po/2 = 0.08 fm"-', in -■"*U+^''^U central 
collisions at a center-of-mass energy £(..„,, = 1200 MeV. Evolutions associated to the four initial 
configurations XX, YX, YY , and ZY are plotted in columns (time runs from top to bottom). Con- 
secutive snapshots are separated by 1.125 zs. Adapted from Ref. [187] 



As an example. Fig. 24 shows snapshots of the density in 238jj^238y central 
collisions at Ec.m. = 1200 MeV. Different shape evolutions are observed depending 
on the initial orientations of the actinides. In particular, a collision of the tips {XX 
configuration) leads to a rapid neck formation, but to a faster re-separation of the 
fragments than with the other orientations. The most compact configurations are 
obtained for side on side collisions {YY and ZY configurations). In particular, the 
ZY configuration leads to the longest contact times as it has less Coulomb repulsion 
than the YY orientation. 

These orientations exhibit also differences as far as the mass flow between the 
di-nuclear fragments is concerned. In fact, no net transfer is observed due to sym- 
metry reasons, except when a tip collide with a side {YX configuration). Indeed, in 
this case, no spatial symmetry prevents a net mass transfer to occur between the 
fragments. In fact, TDHF calculations predict that a large amount of nucleons are 
transferred from the tip to the side, allowing for the production of neutron-rich frag- 
ments in the fermium (Z ~ 100) region in the 238^^238^ j-gaction at energies around 
the Coulomb barrier [52]. This phenomenon has also been investigated with the 
^^^Th-H^^^Cf system using the TDHF approach. In particular, it is shown that when 
the tip of the ^^^Th collide with a side of the ^^"Cf, the latter increases its mass, pro- 
ducing new neutron-rich transfermium nuclei [53, 187]. This phenomenon is called 
"inverse quasi-fission" as the exit channel is more mass asymmetric than the collid- 
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ing partners. Note that inverse quasi-fission is also expected from shell effects in the 
208p|^ region, as shown by calculations based on the Langevin equation [219] (See 
chapter? of Clusters in nuclei Vol. 1 [10]). 



Fig. 25 Collision times for 
each orientation as a function 
of the center-of-mass energy. 
Here, these times are defined 
as the time during which 
the neck density exceeds 
po/10 = 0.016 fm-3. The 
shaded area indicates the 
limit of 2 zs above which 
vacuum decay is expected 
to be observable in central 
collisions. Adapted from 
Ref [52]. 
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Let us now investigate in more details the role of the orientation on collision 
times. Fig. 25 gives the evolution of the contact time between the fragments in 
238y^238y central collisions as a function of energy and for the different orien- 
tations represented in the top of Fig. 24. A saturation of the contact time to 2 zs 
is observed for tip on tip collisions {XX) up to Ec.,,,. ~ 1200 MeV. (Above this en- 
ergy, XX contact times increase due to ternary quasi-fission [52].) This saturation 
is interpreted as a repulsive force generated by large densities in the neck when the 
tips overlap. A similar effect is observed in '^^Th-n^^^Cf [53]. The maximal density 
in the neck region for this system is reported in Fig. 26. We observe that the XX 
configuration leads to large densities above the saturation density. 



Fig. 26 Maximal density in 
the neck as a function of 
the center-of-mass energy in 
^^^Th-F-^"Cf central collisions 
with different orientations 
(see text). The dashed line 
represents the saturation 
density at po = 0.16 fm^-*. 
Adapted from Ref [53] 
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These large densities in the neck are also responsible, in some cases, for a re- 
separation of the system in three fragments [52]. This is illustrated in Fig. 27 for 
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Fig. 27 Zeptosecond snapshots of ^■'**U+~^^*U at £[,„i. = 900 MeV for a tip on tip central collision. 
Adapted from Ref. [52]. 

the 238y^238y Qentral tip on tip collision at Ecm. = 900 MeV. A small fragment is 
observed at rest in the exit channel. A similar phenomenon is discussed in Ref. [10]. 
The formation of this third fragment is interpreted as due to an excess of density 
in the neck region. To illustrate this argument. Fig. 28 shows the internal density at 
the distance of closest approach. Densities above the saturation density are indeed 
observed in the neck region. In fact, instead of breaking in the middle of the neck, 
the system breaks in both sides of this over saturation density region, producing a 
third small fragment at rest [52]. 

Fig. 28 Nucleon density (in 
fm^^) in the collision plane 
is plotted when the density in 
the neck reaches its maximum 
in the XX configuration at 
Ecm. = 900 MeV. The half 
cut surface is an isodensity at 
half the saturation density, i.e. 
Po/2 = 0.08 fm-^. Adapted 
from Ref. [52]. 
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To conclude, the dynamics of actinide collisions has been studied with the TDHF 
approach. These systems encounter quasi-fission with typical di-nuclear system life- 
times of 2 to 4 zs, depending on the initial orientation of the nuclei. The mass 
transfer also strongly depends on the relative orientation at contact. In particular, 
collisions where the tip of a nucleus is in contact with the side of its collision part- 
ner lead to important mass transfer. These multi-nucleon transfer reactions may be 
used in the future for the production of neutron-rich transfermium nuclei. Finally, 
it is shown that the complex neck dynamics may lead to the production of a third 
fragment at rest. 



8.2 Lifetime and spontaneous e^e emission 

It is predicted that nuclear systems with more than ^ 173 protons generate super- 
critical Coulomb fields [10, 220, 221, 222], i.e., with the lowest quasi-molecular 
electronic state in the Dirac sea. This is illustrated in Fig. 29 in the case of a collision 
of two bare uranium. Without Coulomb field, the lowest energy of an electron is 
E = nigC^ for an electron at rest. In the top of Fig. 29, we see that when the two 
uranium are at a distance Du^u ~ 300 fm, the lowest quasi-molecular state available 
for an electron has an energy £ ~ [220]. At this distance, the system is still sub- 
critical as the state is above the Dirac sea. 

A super-critical state is obtained at Du-u — 36 fm when the quasi-molecular 
state crosses the Dirac sea at £ = —meC^ (middle panel in Fig. 29). At this distance 
the nuclei are not yet in contact. The latter occurs at Z)f/-j/ ~ 16 fm. At this distance, 
the quasi-molecular state has an energy £2^—1 .SnieC^. 

If the super-critical state is not or partially occupied, then it induces a hole in 
the Dirac sea. According to QED, such a state is unstable, i.e., it is a resonance 
with a finite lifetime depending on the depth of the energy level. This resonance is 
predicted to decay by producing a e+e^ pair (see bottom panel in Fig. 29). In this 
case, the occupies the state with negative energy E, while the positron is emitted 
in the continuum with an energy —E. 

The Dirac equation has been solved to predict the lifetime of the resonance 
in Ref. [222]. Fig. 30 shows the resulting positron spectra for different hypothe- 
ses about the contact time T between the actinides. The positrons emitted in the 
case r = (solid red line) are due to standard e+e^ pair creation when electric 
charges are accelerated in a strong Coulomb field. For finite contact times, additional 
positrons are produced from the decay of the vacuum. These additional positrons 
form a peak around £ = 1.5 — 2mc^ which can be observed for contact times greater 
than 2 zs. 

To enable the observation of the QED vacuum decay via e+e^ pair emission, it is 
then crucial to have contact times between the actinides as large as possible. Predic- 
tive calculations of the nuclear dynamics, in particular of the contact times between 
the nuclei are then mandatory. We see in Fig. 25 that contact times greater than 2 zs 
are obtained for central collisions at £c.m. > 1000 MeV. In particular, the largest con- 
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Fig. 29 Coulomb potential energy (blue solid lines) generated by two fully stripped uranium nuclei 
at distances Du-u = 299, 31 and 16 fm. The horizontal purple and red lines indicate the lowest 
quasi-molecular electronic state [220]. 
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Fig. 30 Total differential 
positron spectrum dP/clE 
for the U-U system at £(■.,„, = 
740 MeV, for different nuclear 
sticking times calculated 
using the time-dependent 
Dirac equation. The solid red 
curve is for T = zs, the 
long-dashed green curve for 
r = 1 zs, the short-dashed 
blue curve for T = 2 zs, the 
dotted magenta curve for 
r = 5 zs, and the dash-dotted 
cyan curve for T = 10 zs. 
Adapted from Ref [222]. 

tact times of 3-4 zs obtained at Ec.,,,. ~ 1200 MeV might enable the observation of 
spontaneous QED vacuum decay. Similar conclusions have been reached with quan- 
tum molecular dynamics calculations [225] . Note that these energies are greater than 
the previous experimental search for spontaneous e+e^ emissions. Indeed, these 
experiments were using beams around the Coulomb barrier at 700-800 MeV (see 
Ref. [227] and references therein). 

To sum up, a possible signature of the QED vacuum decay would be the obser- 
vation of a peak in the spectra of positrons produced in actinide collisions. Micro- 
scopic calculations predict optimum center-of-mass energies at more than 1 GeV 
These energies are much larger than the energies considered in early experiments. 



9 Dynamics of neutron star crust 

Although the description of supernova explosion mechanisms is not yet complete 
[228], it is well known that "neutron stars" are a possible residue of such explo- 
sions. These dense baryonic objects have drawn lots of interests in the past [229]. 
However, their structure, which is intimately linked to the nuclear phase diagram, 
is still actively discussed nowadays. For instance, the order of the crust-core phase 
transition is under debate [230, 231]. In addition, the structure of the inner core may 
be affected by possible phase transitions involving strangeness [232]. 

The outer layers of neutron stars are also expected to exhibit exotic structures 
involving different geometrical organisations of the nuclear species [234]. Fig. 31 
provides a schematic picture of neutron stars [233]. The outer crust is expected to 
exhibit a lattice of droplets in a sea of nuclear matter [235, 236, 237]. In this picture, 
the inner crust is composed of nuclear clusters which may be unstable to quadrupole 
deformations, forming prolate nuclei. The latter would eventually join up to form 
stringlike structures in the so-called "nuclear pasta phase" [238]. 

The total energy is shared between bulk, surface, and Coulomb contributions. 
Inside nuclei, the density is approximatively constant. As a result, only the surface 
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Fig. 31 Schematic picture of the structure of a neutron star as a function of density. From 
Ref. [233] 



and Coulomb terms depend on nuclear shapes. The details of the diagram in Fig. 31 
are then sensitive to a subtle competition between the surface tension, which makes 
nuclei spherical, versus the Coulomb force, which tends to deform them. 



Fig. 32 Neutron density 
profiles at a temperature 
T = 2.5 MeV, a proton frac- 
tion Xp = 0.3, and baryon 
densities of 0.04 fm^-* (top 
left), 0.06 fm-3 (top mid- 
dle), 0.08 fm-^^ (top right), 
0.09 fm-3 (bottom left), 
0.10 fm-3 (bottom middle), 
and 0.1 1 fm"^ (bottom right) 
The lowest neutron densi- 
ties are in dark (blue) colors, 
while the highest are in gray 
(red). From Ref. [239] 

Static Hartree-Fock calculations based on Skyrme functionals have been per- 
formed to investigate these structures in more details in Refs. [240, 241, 239]. As 
an example. Fig. 32 represents neutron densities for different total baryon densities 
from 0.04 to 0.1 1 fm^^ computed by Newton and Stone [239]. Except at the highest 
density, where an homogeneous distribution is obtained, various pasta phases can 
be observed. 
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The formation of these structures has been also investigated within microscopic 
theories [242, 243, 62, 63]. For instance, the quantum molecular dynamics (QMD) 
semi-classical model has been used by Watanabe and collaborators to show that 
pasta phase could be formed dynamically [242, 243]. 



Fig. 33 Snapshots of the 
density profiles of a per- 
turbed supercell of oxygen 
isotopes in a face-centered cu- 
bic lattice. The mean neutron 
density is (p) = 0.058 fm-^ 
and a proton fraction Xp = 0.5 
is considered. From Ref. [63] 
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In addition, using a quantum microscopic framework based on the TDHF equa- 
tion (called the DYWAN model), Sebille and collaborators have investigated the 
stability of some lattice structures [62, 63]. As an example, the evolution of an ini- 
tially perturbed face-centered oxygen lattice with proton fraction .x:^ =0.5 and mean 
neutron density (p) = 0.058 fm^^ toward a cylinder is shown in Fig. 33. The initial 
perturbation consists of a random displacement of the oxygen centers in order to 
break the mean-field symmetry. 



Fig. 34 Neutron structure 
distributions as a function 
of the threshold density p, 
at which isodensities are 
plotted and of the neutron 
mean density (p) normalised 
to Poo = 0.145 fm"^^ for a 
proton fraction Xp = 0.2. 
Initial simple cubic lattices 
of oxygen isotopes are used. 
From Ref. [63] 
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The role of proton fraction, lattice geometry, and nuclear species on the lattice 
sites have also been considered in Ref. [63]. Different shapes of the pasta phase 
have then been obtained, depending on the mean neutron density (p) and on the 
threshold density used to represent the isodensities. An example of distributions 
of structures as a function of these quantities is shown in Fig. 34. We see that the 
TDHF mean-field dynamics is able to reproduce the different types of structures 
expected in pasta phases. 

To conclude, the dynamics of nuclear cluster aggregation and pasta phase for- 
mation from regular lattice of nuclei can be studied at the TDHF level. It should be 
noted, however, that these dynamics may be sensitive to the energy-density func- 
tional, and, in particular, to its isospin and density dependences [244]. In addition, 
the neutron structure distributions could be different if non-local effective forces 
were used [245, 246]. Finally, the role of beyond mean-field correlations remains to 
be investigated. 



10 Selected conclusions and perspectives 

The TDHF approach provides a mean-field description of nuclear dynamics in the 
presence of some clustering and molecular effects. The same formalism is used 
to investigate light systems exhibiting a-clustering, heavy-ion collisions in a wide 
energy and mass range, and instabilities of lattice of nuclei in neutron star crusts. 

Non-trivial dynamics are observed depending on the initial conditions. For in- 
stance, a-clusters, which are shown to survive and vibrate several zeptoseconds in 
light systems such as '*HeH-^Be— !>^^C, quickly dissolve when entering the mean-field 
of a heavy nucleus. Another example is the formation of di-nuclear systems in col- 
lisions of atomic nuclei. These systems are possible doorways to the formation of 
molecular states of the compound system. 

The dynamics of di-nuclear systems is crucial as it determines the outcome of 
the reaction, i.e., a re-separation of the fragments or their fusion in a compound nu- 
cleus. Although they exhibit some common features across the nuclear chart, such 
as lifetimes of the order of few zeptoseconds, the evolution of di-nuclear systems 
may strongly depend on the entrance channel properties. For instance, fusion prob- 
abilities are much larger in light systems, while the heaviest ones formed in actinide 
collisions always encounter quasi-fission. It is also shown that deformation and ori- 
entation are playing an important role in the dynamics. 

In terms of theoretical description of the proton-neutron composition of the final 
fragments, the TDHF approach can be safely used for the less violent collisions, 
such as quasi-elastic transfer reactions, although one should keep in mind that only 
sequential transfer is included, i.e., transfer of clusters of nucleons is neglected. 
However, damped reactions, such as deep-inelastic collisions and quasi-fission, usu- 
ally involve a large transfer of nucleons between the fragments which may be under- 
estimated at the mean-field level. In fact, experimental charge and mass distributions 
in damped collisions are usually wider than TDHF predictions. This drawback is 
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overcame thanks to the inclusion of fluctuations at the TDRPA level. In addition, 
not only particle number fluctuations are important in damped collisions, but also 
correlations between proton and neutron distributions are shown to be large. These 
correlations, which can be computed within the TDRPA approach, are crucial for a 
good description of fragment distributions in damped collisions. 

The transfer of correlated nucleons (paired nucleons, a-clusters...) is one of the 
main challenges for microscopic approaches. For instance, the inclusion of pairing 
correlations at the TDHFB level will help to describe pair transfer reactions. In par- 
ticular, a possible Josephson effect similar to what is observed in superconductors 
could be observed. 

A good understanding of radiative capture between light nuclei at deep sub- 
barrier energies is crucial to describe stellar nucleosynthesis. However, the present 
time-dependent mean-field approaches do not enable a tunneling of the many-body 
wave-function. Beyond TDHF approaches should then be considered to describe 
sub-barrier fusions of astrophysical interest. In particular, dynamical long-range 
fluctuations must be included. Possible candidates are the time-dependent GCM, 
where the generator coordinate is the distance between the fragments [4], and the 
imaginary time-dependent Hartree-Fock formalism derived from the Feynman path- 
integral approach for many -body systems [67]. 
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